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We describe a simple nearest-neighbor Ising model that is capable of supporting a gas, liquid, 
crystal, in characteristic relationship to each other. As the parameters of the model are varied one 
obtains characteristic patterns of phase behavior reminiscent of continuum systems where the range 
of the interaction is varied. The model also possesses dynamical arrest, and although we have not 
studied it in detail these 'transitions' appear to have a reasonable relationship to the phases and 
their transitions. 



In many systems, one can observe a rich interplay be- 
tween phase-separation, critical phenomena, slowed dy- 
namics and solidification, the latter manifesting itself as 
both crystallization and glassification. These questions 
have a long and venerable past, but in modern condensed 
matter theory the range of systems where the issues have 
become relevant 0[Hy] is remarkable. Recent interest in 
model systems extend from colloidal glasses 0, particle 
gels 0, 0, polymeric gels 0, globular protein crystal- 
lization jjjla] and gellation. For example, it has recently 
transpired that even the simplest system with repulsive 
core and short ranged attraction exhibits a large range 
of phase-transitions, dynamical arrestfpossibly with new 
dynamical logarithmic singularity and kinetic phe- 

nomena [ldj . the pertinent control parameter being the 
range of the attraction. So far it has not been possible 
to understand the inter-relation of all these effects. Such 
issues also lie at the heart of many important modern 
problems of materials science, including the formation of 
arrays of particles on optical wavelengths and knowledge- 
based materials design 0, 0] . 

To fully explore these questions we would need to de- 
scribe an extended range of density, across a number of 
different phase boundaries, dealing naturally with criti- 
cality 13], metastabiliy, and arrest phenomena 0, all 
in a coherent fashion, with tools that were applicable and 
reliable across these regimes. In the interim progress can 
be made in various aspects of the problem |8|, llll ll 51 ] . 

The idea that the gas, liquid, crystal, and transitions 
between them can be studied within the same lattice 
model has been raised over the years in some interesting 
studies 0, and complex models have also been stud- 
ied using such simple models [TtI ]. However, there has 
never been any simple (lattice-based) general mechanism 
producing a liquid (considered as a large collections of 
attraction-dominated states with nearly-degenerate en- 
ergies), in appropriate relation to gas and solid. Here, 
beg inning from ideas introduced by Biroli and Mezard 
|l8l Hf3 ] we show for the first time that it is possible to 
caricature all the states, and their transitions, in a re- 
markably simple nearest-neighbor Ising model. In this 
model dynamical arrest and glassy states are also natu- 
rally incorporated into the story. 



In the model, space is divided into cubes of side a, char- 
acteristic of the particle size, and the microscopic length 
of the system. To the center of each cube we associate 
a site (coordination number c = 6, in a tridimensional 
space) and Ising-like occupancy variable. Where a site 
is occupied we define (respectively) attractive and repul- 
sive interactions between c a and c — c r of the nearest 
neighbors. 

Thus c r specifies the fraction of space within a given 
cube that is available to neighboring particles, given the 
nature of the repulsive interactions, whilst c a specifics 
how many of those particles can benefit from the attrac- 
tive energy of the system. Whilst we emphasize that 
this is a caricature, rather than a representation, of the 
potential, we expect that c r therefore expresses the com- 
plexity or irregularity of the core of the particle, reflected 
ultimately in the density at random close packing in the 
real system, and in the model. The interpretation of c a 
is more subtle. However, broadly speaking, a small num- 
ber of (c a ) attractive interactions are shared between a 
larger number (c r ) neighbors, there are many local con- 
figurations in which the attractive contribution saturates 
in the dense system. This freedom encourages the forma- 
tion of a liquid. This idea is consistent with the rationale 
for formation of a gel by short-ranged colloidal particles 
0. In that case there is such limited freedom by which 
the attractions can be captured that the liquid arrests to 
form a gel-like state. There is also a general trend for 
the model to be dominated by the 'liquid' state as we 
approach c r = c a = 6, the nearest-neighbor model. 

In our present study, we assume nearest-neighbor finite 
attractions (e) and infinite repulsions, the latter to make 
contact with related studies |18lll^ |. 

Thus, the infinite repulsion between c r neighbors may 
be written (0 < c r < 6): 
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where lj is the number of occupied nearest-neighbor sites. 
In the case of c r = 6 no repulsive interaction is present in 
the model. The attractive interaction is governed by the 
parameter c a (0 < c a < 6) and is defined (per particle) 
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where e > is the strength of the attraction. Throughout 
the paper we will define energies and chemical potential 
in units of kf,T, and the 'effective' energy will refer to the 
combined exponent of energy and chemical potential in 
the partition function. 

Despite the fact that this is a simple nearest-neighbor 
Ising model, it is remarkably rich and spatial frustra- 
tion produces a large range of phenomena, some of which 
have previously been observed in more complicated lat- 
tice models. In such cases it has also been possible to 
formulate methods for their study HI HHH H3 • For 
the present case we show that the existence of a liquid 
disposed between the gas and liquid, and resulting triple 
point originate in a large degeneracy of states that are 
found at zero temperature as 'multiphase points'. In spe- 
cial cases we make the number of such states macroscopic, 
but even where they are not so, at finite temperatures 
we can cause the existence of an infinite number of near- 
degenerate states, and a liquid. 

It is possible to determine the phase-diagram of such 
models to a high degree of confidence. To begin with 
we determine the ordered phases by construction of the 
zero temperature states of the model. On the simple 
cubic lattice we re-write the underlying effective Hamil- 
tonian (itself a combination of pure energy and chemical 
potential) as a sum over octahedral fragments [U 
to each of which we can associate an independent en- 
ergy. We may then choose the 'optimal' fragments that 
are able to tile the lattice, and thereby classify the 
zero-temperature states. There is at least one multi- 
phase point [2H E2I Eiil at which an infinity of zero- 
temperature states is degenerate, though the entropy per 
particle remains finite. 

In general terms the zero temperature states are or- 
ganized as follows. In three dimensions, there are 20 
different types of fragment. Now, in the Grand Canon- 
ical ensemble, for fixed values of c r and c , the energy 
of a particle-centered fragment has the form — qe — p 1 
where p is the chemical potential and q is a positive inte- 
ger number depending on the number of neighbors. All 
the vacancy-centered fragments have zero energy, unlike 
previous frustrated models. Therefore we use a slightly 
different notation, defining as n and a the total number 
of nearest neighbors and the number of filled axes (lin- 
ear triplets of sites, the extreme pair being occupied), 
respectively, and using a label p = 0, 1 on the left of the 
symbol of a fragment distinguishing the vacancy-centered 
and the particle-centered ones: p f£. 

The internal energy per volume of a state of the model 



where L 3 is the lattice volume and pj is the density of 
the j—th fragment; the sum is intended only over the 
particle-centered fragment types (PC). 

For p — > — 00 the ground state is the zero-density gas, 
because in that limit the energy of every particle-centered 
fragment is infinite. 

For p — » +oo the ground state is the most dense crys- 
tal, i.e. the crystal with the lowest density of vacancy- 
centered fragments. In between we have many different 
situations, depending on the choice of c a and c r . Below 
we give only illustrative examples. 

Thus, for the case c a = 3 and c r = 3 we have two 
fragments with energy — 3e — p, but only one d/3) tiles 
the space with the highest density vacancy-centered crys- 
tal /g. So there is one highest density crystal which is 
ordered along the diagonal with a repeat pattern of an 
occupied pair of layers and an unoccupied layer. This 
state is degenerate with the zero-density gas at the point 
p = — 3e, the zero-temperature boundary between gas 
and the crystal. There cannot be finite regions occupied 
by other states at zero temperature. In fact, a different 
state could have lower energy only if 
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and this is impossible because qj < 3 for every j. 

The point p = — 3e is a multiphase point, because 
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FIG. 1: Pure states for the case c r = 3,4, c a = 3. The 
first is made up of f% fragments (top). There are seven pure 
phases made from the /I fragments. Each of the seven pairs 
constitues two layers of a 4 x 4 x 4 unit cell the remainder 
being the particle-vacancy symmetric ones. 
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two particle-centered fragments df°, 1 fl) and every 
vacancy-centered fragment is degenerate. Some examples 
of crystals at the multiphase point are given in Figure ^ 
For c a — 4 and c r = 3 there are two crystals with the 
highest energy 0.75, made up of the couples (1/4,0/!)' 
(i/l) o/l)) which are represented in Figure |2(a)| and 
|2(b)| respectively. The point /i = — 3e is still a mul- 
tiphase point, but now four particle-centered fragments 
(1/311/3; 1/4: 1/4) are degenerate. For this reason, here 
the multiphase point is much more degenerate. 
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FIG. 4: Monte Carlo Grand ensemble simulation for the case 
c r — 3, c a = 3. Graph of the temperature (/3e) _1 versus 
chemical potential (fife). Upper inset: density vs chemical 
potential graph for k^T — 0.5e. Lower inset: long time be- 
havior for the point at [i = — 2.4e. 
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FIG. 2: Unit cells of the crystals at T — 0. The black spheres 
represent an occupied position on the lattice and the white 
ones represent a vacancy. In both cases the density is p — 
0.75. It transpires that the crystal (a) is more stable (at 
low finite temperature) than (b), and therefore only (a) is 
observed by cooling in the simulation. 
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FIG. 3: Monte Carlo Grand ensemble simulation for the case 
Cr = 4, c a — 3. Graph of the temperature (/3e) _1 versus 
chemical potential (p/e). The (*) are the results of the sim- 
ulation lowering the temperature. The (o) are the results of 
the simulation increasing the chemical potential. The criti- 
cal point is labelled C. The inset shows an isothermal scan of 
density and heat capacity for kbT = 0.6e. 



Besides determining the low temperature limit of the 
models, such considerations allow us to estimate, and 
thereby control, as a function of c a and c r , the number 
of low temperature states where attraction is dominant. 
This allows us ultimately to fix the triple point of the 
system. 

We present examples of Grand Canonical simulations 
(L = 12, < 10 6 MCS). (We remark that we have as yet 
attempted no systematic exploration of dynamical arrest 
except the expected transient arrest phenomena of the 
underlying model (lij)- 

Phase diagrams for c,. = 4, c a = 3 and c r — 3, c a = 3 
are shown in Figure |31 and 0] respectively. 

The first is representative of a conventional gas-liquid- 
crystal phase diagram, and the usual range of attractive 
interactions. Thus the simulations readily reproduce the 
zero temperature states of gas and crystal, and there is 
a significant range where the liquid state is present. In 
the inset to Figure we show the density and heat ca- 
pacity in an isothermal cut of the phase-diagram, just 
above the gas-liquid critical point. Such isotherms, as 
the temperature decreases from this point, exhibit an in- 
creasing peak height of the heat capacity, rising towards 
the critical point. The liquid-crystal phase boundary is 
a conventional first-order phase-transition, that asymp- 
totes towards the fluid-repulsive crystal transition, as ex- 
pected. Simulations near the (low but finite T) triple 
point become quite slow. 

The second phase diagram is chosen to show how the 
model represents a short range potential. The simula- 
tions are still for quite small systems. However, the 
gas-liquid branch, if it exists, is so short to be almost 
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unidentifiable. This would be the case, for example, of a 
continuum square-well potential when the well width is 
less the ten percent of the core size. 

Furthermore, it is intriguing to note the numerous 
other indications of a complex equilibrium, arrest, and ki- 
netics picture, as expected for the equivalent continuum 
picture. In particular, there does appear to be a dense 
fluid-like state that is long lived, as well as a string-like 
low density fluid that may credibly identified as a gel, 
both expected. (In the upper inset to the figure we show 
an isothermal scan in which the system makes an ap- 
parent transition to a dense liquid-like phase, whereas in 
the lower inset we show that this state degenerates to a 
crystal at very long times) 0, 0] ■ 

In summary, we have shown that a very simple Ising 
lattice model, with only nearest neighbor interactions, 
is capable of supporting all of the characteristic phenom- 
ena of particles with repulsive and attractive interactions; 
gas, liquid, crystals, 'glass-transition' phenomena of var- 
ious kinds, and that variation of parameters provides a 
mechanism to mimic the range of the potential into im- 
portant regimes that have yet to be understood. 

Many of the outstanding problems of interest such as 
metastability or buried critical points, arrest phenomena, 
and near-criticality and gellation are expected to be quite 
well-described on lattices. Crystallization nucleii and to 
some degree their growth are the most dubious aspects of 
a lattice description, but even there many of the most dif- 
ficult and subtle unexplained processes currently involve 
some element of near-criticality |l3j . something which is 
well-dealt with in the lattice, and is difficult to deal with 
in a useful manner otherwise. 

Besides having a certain elegance in being a simplest 
model known to exhibit all of the relevant phenomena, 
it seems likely that this model has considerable relevance 
to important practical issues. 
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